Evidence for retrograde lithospheric subduction on Venus

Though there is no plate tectonics per se on Venus, recent Magellan radar images and topographic profiles of the planet suggest the occurrence of the plate tectonic processes of lithospheric subduction and back-arc spreading. The perimeters of several large coronae (e.g., Latona, Artemis, and Eithinoha) resemble Earth subduction zones in both their planform and topographic profile. The planform of arcuate structures in Eastern Aphrodite were compared with subduction zones of the East Indies. The venusian structures have radii of curvature that are similar to those of terrestrial subduction zones. Moreover, the topography of the venusian ridge/trench structures is highly asymmetric with a ridge on the concave side and a trough on the convex side; Earth subduction zones generally display the same asymmetry

Bounds on Lithospheric Thickness on Venus from Magellan Gravity and Topography Data

The primary objective of the work executed under NAGW-4784 is to provide constraints on the thermal and tectonic evolution of Venus. Establishing thermal and tectonic evolution models requires not only geological, but geophysical constraints, in particular the nature of temporal and spatial variations in crustal and lithospheric thickness. The major topics of study completed under NAGW-4784 (described more fully below) are: (1) detailed analyses of the resolution of Magellan Line-Of-Site (LOS) Doppler data to establish the minimum resolvable wavelength in the gravity data; (2) calculations of the global strain field in the venusian lithosphere and comparisons with global strain patterns from geological mapping; (3) study of the geological history of coronae at E. Eistla Regio; (4) estimation of crustal and lithospheric thickness by modeling of topography at asymmetric and symmetric rift-like chasmata; (5) preliminary investigations of spatial versus temporal variations in lithospheric thickness. Both the PI and Co-I have presented papers based on these topics at national and international meetings (American Geophysical Union Meetings, Lunar and Planetary Science Conferences, Chapman Conference on the Geodynamics of Venus)

Experimental Design for Bathymetry Editing

We describe an application of machine learning to a real-world computer assisted labeling task. Our experimental results expose significant deviations from the IID assumption commonly used in machine learning. These results suggest that the common random split of all data into training and testing can often lead to poor performance.Comment: Published as a workshop paper at ICML 2020 Workshop on Real World Experiment Design and Active Learnin

Outer Trench Slope Flexure and Faulting at Pacific Basin Subduction Zones

Flexure and fracturing of the seafloor on the outer trench wall of subduction zones reflects bending of the lithosphere beyond its elastic limit. To investigate these inelastic processes, we have developed a full non-linear inversion approach for estimating the bending moment, curvature, and outer trench wall fracturing using shipboard bathymetry and satellite altimetry derived gravity data as constraints. Bending moments and downward forces are imposed along curved trench axes and an iterative method is used to calculate the non-linear response for 26 sites in the circum-Pacific region having seafloor age ranging from 15 to 148 Ma. We use standard thermal and yield strength envelope models to develop the non-linear moment versus curvature relationship. Two coefficients of friction of 0.6 and 0.3 are considered and we find the lower value provides a better overall fit to the data. The main result is that the lithosphere is nearly moment saturated at the trench axis. The effective elastic thickness of the plate on the outer trench slope is at least three times smaller than the elastic thickness of the plate before bending at the outer rise, in agreement with previous studies. The average seafloor depth of the unbent plate in these 26 sites matches the Parsons & Sclater (1977) depth versus age model beyond 120 Ma. We also use the model to predict the offsets of normal faults on the outer trench walls and compare this with the horst and graben structures observed by multibeam surveys. The model with the lower coefficient of friction fits the fault offset data close to the trench axis. However, the model predicts significant fracturing of the lithosphere between 75 and 150 kilometres away from the trench axis where no fracturing is observed. To reconcile these observations, we impose a thermoelastic pre-stress in the lithosphere (Wessel 1992) prior to subduction. This pre-stress delays the onset of fracturing in better agreement with the data

Report of the panel on lithospheric structure and evolution, section 3

The panel concluded that NASA can contribute to developing a refined understanding of the compositional, structural, and thermal differences between continental and oceanic lithosphere through a vigorous program in solid Earth science with the following objectives: determine the most fundamental geophysical property of the planet; determine the global gravity field to an accuracy of a few milliGals at wavelengths of 100 km or less; determine the global lithospheric magnetic field to a few nanoTeslas at a wavelength of 100 km; determine how the lithosphere has evolved to its present state via acquiring geologic remote sensing data over all the continents

insar decorrelation to assess and prevent volcanic risk

SAR� can� be� invaluable� describing� pre�eruption� surface� deformation� and� improving� the� understanding� of� volcanic� processes.� This� work� studies� correlation� of� pairs� of� SAR� images� focusing� on� the� inༀ䃻uence� of� surface,� climate� conditions� and� acquisition� band.� Chosen� L�band� and� C�band� images� (ENVISAT,� ERS� and� ALOS)� cover� most� of� the� Yellowstone� caldera� (USA)� over� a� span� of� 4� years,� sampling� all� the� seasons.� Interferograms� and� correlation� maps� are� generated� and� studied� in� relation� to� snow� depth� and� temperature.� To� isolate� temporal� decorrelation� pairs� of� images� with� the� shortest� baseline� are� chosen.� Results� show� good� performance� during� winter,� bad� attitude� towards� wet� snow� and� good� coherence� during� summer� with� L�band� performing� better� over� vegetation

Oceanic microplate formation records the onset of India-Eurasia collision

Mapping of seafloor tectonic fabric in the Indian Ocean, using high-resolution satellite-derived vertical gravity gradient data, reveals an extinct Pacific-style oceanic microplate ('Mammerickx Microplate') west of the Ninetyeast Ridge. It is one of the first Pacific-style microplates to be mapped outside the Pacific basin, suggesting that geophysical conditions during formation probably resembled those that have dominated at eastern Pacific ridges. The microplate formed at the Indian-Antarctic ridge and is bordered by an extinct ridge in the north and pseudofault in the south, whose conjugate is located north of the Kerguelen Plateau. Independent microplate rotation is indicated by asymmetric pseudofaults and rotated abyssal hill fabric, also seen in multibeam data. Magnetic anomaly picks and age estimates calculated from published spreading rates suggest formation during chron 21o (~47.3 Ma). Plate reorganizations can trigger ridge propagation and microplate development, and we propose that Mammerickx Microplate formation is linked with the India-Eurasia collision (initial 'soft' collision). The collision altered the stress regime at the Indian-Antarctic ridge, leading to a change in segmentation and ridge propagation from an establishing transform. Fast Indian-Antarctic spreading that preceded microplate formation, and Kerguelen Plume activity, may have facilitated ridge propagation via the production of thin and weak lithosphere; however both factors had been present for tens of millions of years and are therefore unlikely to have triggered the event. Prior to the collision, the combination of fast spreading and plume activity was responsible for the production of a wide region of undulate seafloor to the north of the extinct ridge and 'W' shaped lineations that record back and forth ridge propagation. Microplate formation provides a precise means of dating the onset of the India-Eurasia collision, and is completely independent of and complementary to timing constraints derived from continental geology or convergence histories. © 2015 Elsevier B.V.K.J.M. and R.D.M. were supported by ARC Discovery Project DP130101946 . The CryoSat-2 data were provided by the European Space Agency, and NASA/CNES provided data from the Jason-1 altimeter. This research was supported by the National Science Foundation ( OCE-1128801 ), the Office of Naval Research ( N00014-12-1-0111 ), the National Geospatial-Intelligence Agency ( HM0177-13-1-0008 ) and ConocoPhillips . Version 23 of the global grids of gravity anomaly and VGG can be downloaded from the following ftp site ftp://topex.ucsd.edu/pub/global_grav_1min . All figures were produced using the Generic Mapping Tools ( GMT ) software ( Wessel et al., 2013 ). The open-source plate reconstruction software GPlates ( Boyden et al., 2011 ) was used to compute the distance from the Kerguelen Plume to the initiation point of ridge propagation using different absolute reference frames, and to produce the VGG raster reconstruction in Fig. 3 . Magnetic anomaly picks were accessed from the compilation of Seton et al. (2014) from The Global Seafloor Fabric and Magnetic Lineation Data Base Project website ( http://www.soest.hawaii.edu/PT/GSFML/ ). We thank the editor An Yin and two anonymous reviewers for their thoughtful and constructive comments that improved the manuscript. Appendix

Refining the shallow slip deficit

Geodetic slip inversions for three major (M_w > 7) strike-slip earthquakes (1992 Landers, 1999 Hector Mine and 2010 El Mayor–Cucapah) show a 15–60 per cent reduction in slip near the surface (depth < 2 km) relative to the slip at deeper depths (4–6 km). This significant difference between surface coseismic slip and slip at depth has been termed the shallow slip deficit (SSD). The large magnitude of this deficit has been an enigma since it cannot be explained by shallow creep during the interseismic period or by triggered slip from nearby earthquakes. One potential explanation for the SSD is that the previous geodetic inversions lack data coverage close to surface rupture such that the shallow portions of the slip models are poorly resolved and generally underestimated. In this study, we improve the static coseismic slip inversion for these three earthquakes, especially at shallow depths, by: (1) including data capturing the near-fault deformation from optical imagery and SAR azimuth offsets; (2) refining the interferometric synthetic aperture radar processing with non-boxcar phase filtering, model-dependent range corrections, more complete phase unwrapping by SNAPHU (Statistical Non-linear Approach for Phase Unwrapping) assuming a maximum discontinuity and an on-fault correlation mask; (3) using more detailed, geologically constrained fault geometries and (4) incorporating additional campaign global positioning system (GPS) data. The refined slip models result in much smaller SSDs of 3–19 per cent. We suspect that the remaining minor SSD for these earthquakes likely reflects a combination of our elastic model's inability to fully account for near-surface deformation, which will render our estimates of shallow slip minima, and potentially small amounts of interseismic fault creep or triggered slip, which could ‘make up’ a small percentages of the coseismic SSD during the interseismic period. Our results indicate that it is imperative that slip inversions include accurate measurements of near-fault surface deformation to reliably constrain spatial patterns of slip during major strike-slip earthquakes

Deformation on Nearby Faults Induced by the 1999 Hector Mine Earthquake

Interferometric Synthetic Aperture Radar observations of surface deformation due to the 1999 Hector Mine earthquake reveal motion on several nearby faults of the eastern California shear zone.We document both vertical and horizontal displacements of several millimeters to several centimeters across kilometerwide zones centered on pre-existing faults.Portions of some faults experienced retrograde (that is, opposite to their long-term geologic slip) motion during or shortly after the earthquake.The observed deformation likely represents elastic response of compliant fault zones to the permanent co-seismic stress changes. The induced fault displacements imply decreases in the effective shear modulus within the kilometer-wide fault zones, indicating that the latter are mechanically distinct from the ambient crustal rocks